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Melting of Single Lipid Components in Binary 
Lipid Mixtures: A Comparison between FTIR 
Spectroscopy, DSC and Monte Carlo Simula- 
tions 



M. Fidorra^, T. Heimburg^ and H.M. Seegert'*'* 

t The Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark 
■f MEMPHYS-Center of Biomembrane Physics, Department of Physics, 
University of Southern Denmark, Odense, Denmark 
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Modena, Italy 

Monte Carlo (MC) Simulations, Differential Scanning Calorimetry (DSC) and Fourier Transform InfraRed (FTIR) 
spectroscopy were used to study the melting behavior of single lipid components in two-component membranes of 1,2- 
Dimyristoyl-D54-sn-Glycero-3-Phosphocholine (DMPC-d54) and l,2-Distearoyl-sn-Glycero-3-Phosphocholine (DSPC). 
Microscopic information on the temperature dependent melting of the single lipid species could be investigated using FTIR. 
The microscopic behavior measured could be well described by the results from the MC simulations. These simulations 
also allowed to calculate heat capacity profiles as determined with DSC. These ones provide macroscopic information about 
melting enthalpies and entropy changes which are not accessible with FTIR. Therefore, the MC simulations allowed us to 
link the two different experimental approaches of FTIR and DSC. 

keywords: domain formation; phase separation; FTIR; MC simulation; DSC; melting of single lipid components 



Introduction 

Biological membranes display a complex lipid composition. 
The reason for the big variety in lipid species is still not clear. 
It has been shown for many different bacteria that their lipid 
synthesis depends on physical parameters such as growth tem- 
perature or hydrostatic pressure (llj,|2j,|^|j,|5j|3|7l|8|,|9l LLOl LLil 
regard also the citations given therein) . It could be shown 
that the membranes of bacteria grown at higher temperatures 
contain more saturated fatty acid chains and chains tend to be 
longer fflllllllllllBIIGilsee also citations therein). 
Higher hydrostatic pressure induces an increased synthesis of 
unsaturated lipids (|9l llOl) . which has also been found for bi- 



ological membranes from several deep-sea fish tissues (1121) . 
This adaption has an influence on the physical behavior of a 
biological membrane such as the melting transition behavior. 
Melting transitions from a solid ordered (so; in the litera- 
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sn-glycero-3-phosphocholine 
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ture also called gel) to a liquid disordered {Id; or named as 
fluid) phase are present in artificial and biological membranes 

(E til H EH E3, H). In the presence of cholesterol also 

liquid ordered phases develop (1191) . Considering single lipid 
membranes one finds that a higher degree of saturation or 
longer fatty acids result in higher melting temperatures J20h . 
Therefore, one would expect a change of the melting behav- 
ior of biological membranes in dependence of growth tem- 
perature or hydrostatic pressure and therewith in dependence 
of lipid synthesis. Indeed, this has been demonstrated for 
the prominent example of the inner membrane of Escherichia 
coUdH. 

The dependence of lipid synthesis by outer physical param- 
eters such as temperature or hydrostatic pressure was termed 
"homeoviscous adaption" (I22i 1231) . The underlying idea was 
that for a proper cell function the membrane needs to provide 
a certain viscosity. In this context the widely used term "flu- 
idity" was introduced in the literature. It should, however, be 
noted that this term lacks a clear definition. It was also ar- 
gued against a role of "fluidity" (1241) . An alternative role of 
the adaption of lipid synthesis on outer physical parameters 
is to ensure a control of the heterogeneity of the lateral mem- 
brane structure (25). This idea has obtained special attention 
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within the discussion about "rafts" (26). 

Artificial lipid membranes consisting of one, two or three 
lipid components are suitable systems to study the general 
physical behavior of lipid phase transitions in more complex 
lipid systems . In the phase coexistence regime domains with 
a size ranging from a few nanometers to several micrometers 
form in model systems and can be observed by a large vari- 
ety of experimental techniques, such as Fluorescence Reso- 
nance Energy Transfer, Atomic Force Microscopy and Con- 
focal Fluorescence Microscopy (|27, 28, 2§l 30l [3lb . 

Already since the early 1970s the idea that domain forma- 
tion processes can influence protein activity and function was 
expressed ^251 l32l [33b- In several publications the existence 
of domains have been shown to provide a mechanism to con- 
trol biochemical reaction cascades d34l l35l l36l l37l [38b . Ev- 
idence that lipid phase transitions and domain coexistence 
trigger enzyme activity has been found for various systems. 
For example, the enzyme phospholipase A2 which hydrolyzes 
the sn-2 bond of phospholipids, producing a free fatty acid 
and a lysolipid, is strongest in activity when the lipid mem- 
brane is in the lipid phase transition regime (39, 40, 4l|). The 
activity of the special class of phospholipase A2 type II A has 
been demonstrated to be related to an enrichment of anionic 
lipids into Id domains d42l) . In this case, a different partition 
behavior of the anionic lipids into Id and lo domains led to 
a local concentration higher than the threshold value for en- 
zyme activity. Changes in the opening times and probabilities 
of calcium channels reconstituted in POPE:POPC membrane 
of various ratios have been interpreted to be due to the coex- 
istence of solid ordered and liquid disordered phases d43l) . A 
further example is the protein kinase C (PKC) which modifies 
proteins by chemically adding phosphate groups. It has been 
discussed that lateral heterogeneities of the lipid membrane 
control the activation of this enzyme d44l I45D . In the latter 
study it has also been pointed out that domains enriched in 
dioleoylglycerol lipids play a crucial role, but also a linear re- 
lationship between PKC activity and phosphatidylserine con- 
tent has been found (46). Further communications can be re- 
trieved where a connection between enzyme activity and the 
presence of specific lipids and structure are made. This is e.g. 
also the case for the calcium ATPase ( 147ft . 

The publications cited above are only a small fraction of 
many in that respect, but they already show that it is impor- 
tant not only to study domain formation, but also to investi- 
gate the fraction of the single lipid species in lipid membranes 
showing lateral inhomogeneity due to phase separation. The 
single lipid species will occur in different amounts in the dif- 
ferent domains that form in a phase transition regime (1481) . 
The work presented here concentrates on this for the case of 
a binary lipid mixture composed of DMPC and DSPC. These 
two phospholipids display melting transition events with tran- 
sition midpoint temperatures of 23.6°C or 54.4°C, respec- 
tively (dot). In our study we used a deuterated DMPC derivate 
whose phase transition midpoint temperature lies at around 



19 - 2o o c (SHU. 

For the description of phase separation in lipid membranes 
various theoretical approaches have been applied. This in- 
cludes ideal and regular solution theory, mean field, but also 
numerical approaches such as Monte Carlo (MC) simulations 
d48l pll l52l l53l p4h . Ideal and regular solution theories al- 
low the understanding of phase diagrams. In the ideal solu- 
tion theory complete miscibility in the solid ordered and liq- 
uid disordered phase is assumed. The regular solution theory 
takes a step further and regards a non-ideal mixing in the solid 
ordered phase. These considerations allow the construction 
of phase diagrams. Phase diagrams are used to map the phase 
behavior of lipid membranes in dependence of e.g. temper- 
ature and lipid composition. The theoretically derived lever 
rule provides the possibility to determine the ratio of disor- 
dered and ordered lipids in dependence of temperature and 
lipid composition of binary lipid mixtures and can also be ap- 
plied to experimentally measured phase diagrams (48). Ideal 
or regular solution theory are, however, based on simplify- 
ing assumptions. This does not only include an ideal mixing 
in the liquid phase, but also the assumption of phase separa- 
tion under all conditions in the melting regime. This is, how- 
ever, not the case (1551) . Further, fluctuations which are en- 
hanced in the melting transition regime are present (13, 57). 
This is where MC simulations analyzing simple models of 
lipid chain melting enter. They allow the consideration of in- 
teractions between different lipid species and varying chain 
state. The formation of lipid domains and phases and the 
consideration of fluctuations is a direct consequence of the 
thermodynamics of these systems which is included in these 
models. Therefore, these simulations can describe lipid mem- 
brane properties in the vicinity of lipid phase transitions, such 
as changes in membrane permeabilit y d58h. fluctuations in en- 
thalpy (1571) and diffusion processes (59, 60) well. 

In this study we determined the amount of the single lipid 
species in the different domains of binary DMPC and DSPC 
mixtures by means of MC simulations and compared the re- 
sults obtained to experimental data measured by Fourier Trans- 
form InfraRed spectroscopy (FTIR). 

Among other techniques FTIR is well suited to investigate 
the melting behavior of lipid membranes. In these studies the 
temperature dependence of various vibrational modes due to 
structural changes provide a versatile tool in following phase 
transitions Articles have been published on experiments 
in which FTIR was used to measure nanoscale domain size 
d62D and chain order parameters d63ft . 

In our study the hydrogen atoms of DMPC were replaced by 
deuterium atoms. Due to this the melting of the DMPC lipid 
chains could be separated from the DSPC lipid chain melt- 
ing, as the CH2 and CD2 symmetric and asymmetric stretch 
vibrations differ in frequency due to the heavier deuterium 
atoms. This allowed us to compare the melting of single com- 
ponents with MC simulation results which contain informa- 
tion of the disordered chain ratio of single lipid components. 
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MC simulations are able to describe experimental heat capac- 
ity profiles well. Therefore, they provide a means to couple 
the two different experimental approaches. 

The experimental absorption spectra were evaluated using 
two methods. One of the methods probed the temperature de- 
pendence of the position of the absorption maxima of the anti- 
symmetric CD2 and symmetric CH2 stretch vibrations. The 
other one evaluated peak area changes of difference spectra 
(Difference Spectra method). Simulation results on the frac- 
tion of disordered chains of the single lipid species agreed 
well with the data found using the latter method. FTIR al- 
lows to distinguish the melting of the single lipid species and 
it probes the microscopic behavior. DSC probes the bulk be- 
havior and macroscopic observables. The Monte Carlo simu- 
lations were used as an approach to combine the informations 
of the two experimental techniques. DSC, MC and FTIR are 
a potent combination to investigate the melting behavior of 
lipid membranes. 

Materials and Methods 

Sample Preparation: 

The lipids l,2-dimyristoyl-d54-sn-glycero-3-phosphocholine 
(DMPC-d54) and l,2-distearoyl-sn-glycero-3-phosphocholine 
(DSPC) were purchased from Avanti Polar Lipids (Alabaster, 
USA). Samples were prepared from lipid stock solution in 
chloroform. These were prepared in chloroform from lipid 
powder as delivered from the vendor without further purifi- 
cation. To prepare lipid samples in buffer, chloroform from 
the desired amount of lipid stock solution mixture was evap- 
orated under a stream of nitrogen and the sample was kept 
under vacuum for several hours afterwards. After this pro- 
cedure the dried lipid film was hydrated with preheated ul- 
trapure water at a temperature of 60°C, which is higher than 
the melting temperature of the higher melting lipid compo- 
nent DSPC (54.4° C, (H^)) and kept stirring at this tempera- 
ture for about 30 minutes. In calorimetric experiments a lipid 
concentration of 10 mM was used. The Fourier Transform 
Infrared spectroscopy measurements required concentrations 
of 65 mM. In order to facilitate lipid membrane hydration 
under these high concentrations the samples were subjected 
to several freeze-thaw cycles. 

DSC experiments: 

Differential Scanning Calorimetry (DSC) experiments were 
performed on a VP-DSC (MicroCal, Northampton, MA, USA) 
with a scan rate of 15°C/h (DMPC-d54:DSPC 50 : 50) or of 
5°C/h (DMPC-d54) per hour in upscan direction. The sam- 
ples were equilibrated at the starting temperature for 30 min- 
utes. 

Fourier Transform InfraRed spectroscopy experiments: 

Fourier Transform InfraRed (FTIR) spectroscopy experiments 
were performed on a Vertex70 (Bruker Optics, Bremen, Ger- 



many) equipped with a mercury-cadmium-telluride (MCT) 
detector. The spectrometer was flushed with nitrogen start- 
ing about 30 minutes before the experiments. Temperature 
control was done by a self built sample holder connected to a 
water bath (DC30-K20, ThermoHaake, Karlsruhe, Germany). 
The scan rate for the temperature ramp was set to 20°C/h. The 
temperature in the sample corresponding to each temperature 
point set with the water bath was appointed in a reference 
scan. Background spectra were acquired at each temperature 
point after lifting the sample compartment out of the beam 
path. We were interested in the IR absorption due to the 
symmetric and asymmetric stretch vibrational modes of the 
CH2 and CD2 groups of the lipid fatty acid chains. Example 
absorption maxima of the CH2 symmetric and asymmetric 
stretch vibrations at different temperatures are displayed in 
the left panel of fig. [2] 

Two methods to analyze the FTIR spectroscopy data were 
used. In a first approach FTIR spectroscopy data was pro- 
cessed with a self- written plugin for the IGOR software pack- 
age (WaveMetrics, Lake Oswego, OR, USA). The absorption 
peaks were fitted with a polynomial function. Then the peak 
position was determined by calculating the root of the fit. 

For the second protocol the peak area of the CH2 or respec- 
tively of the CD2 stretch vibrations was normalized to a value 
of 1 x cm . A reference spectrum was chosen which was 
taken at a temperature below the melting transition regime. 
This spectrum was subtracted from all other spectra which 
were recorded at different temperatures (see the right panel 
of fig. [2j. Then the peak area located at the positive axis of 
ordinates was calculated and plotted as a function of temper- 
ature. There was no difference in surveying both, symmetric 
and antisymmetric modes or the single ones. Further, we shall 
call this method the Difference Spectra method. 

Model: 

We have already applied an Ising-like model to the study 
of diffusion J59I) . fluctuation J57I) and relaxation properties 
(1641) of DMPCDSPC lipid mixtures. In this paper we used 
this model to describe FTIR measurements on mixtures of 
DMPC-d54 and DSPC and linked these with DSC measure- 
ments. 



The model has been presented in detail elsewhere (154 15£ 
In short, the model is based on the assumption that the lipid 
chain states can be described by assigning either an ordered 
or disordered state to the fatty acid lipid chains. Lipid chains 
are arranged on a triangular lattice and only nearest neighbor 
interactions are included. Lipid chains might change the de- 
gree of order or their positions during the simulation. The 
model is evaluated with Monte Carlo (MC) simulations in 
which the free energy difference of an old and a trial con- 
figuration is sufficient for decision making of the acceptance 
of a new configuration. The free energy of the system can be 
divided into a configuration dependent and an independent 
part. The configuration dependent contribution equals (54): 
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AG(S) = Nf(AHt - TASi) + N^(AH 2 - AS 2 ) 

+ Nguft + JVgwg + Ngu$ (1) 

+ + + Ngu° 2 i 

where AHi and ASi are the enthalpy and respective entropy 
differences of the disordered to the ordered state of DMPC 
(i = 1) and DSPC (i = 2) lipid chains. uf^ n are inter- 
action parameters of unlike nearest neighbors and N™ n the 
corresponding number of the unlike nearest neighbor con- 
tacts with (i ^ j) and (m ^ n). m and n denote the state 
of a chain which is either ordered (o) or disordered (d). En- 
thalpy, entropy and nearest neighbor interaction parameters 
can be obtained from experimentally determined heat capac- 
ity profiles. Previously, we have determined the parameters 
needed for the DMPCDSPC system i5% . In this paper, how- 
ever, FTIR experiments were performed with DMPC-d54 and 
DSPC mixtures. In the following it was assumed that the un- 
like species interaction parameters did not change from the 
system DMPCDSPC to DMPC-d54:DSPC. The heat capac- 
ity profile of the pure DMPC-d54 was measured to obtain its 
enthalpy and entropy changes during the melting transition 
and the cooperativity parameter u)°f {data not shown). The 
remaining parameters were taken from a recent publication 
j59l) . All values are listed in table Q] 
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Table 1 : Parameter values used in the Monte Carlo simulations. They were 
determined from experimentally determined heat capacity profiles 1: 591) . All 
numbers are given per lipid chain. The indices o and d stand for ordered and 
disordered, respectively. 1 and 2 index DMPC-d54 and DSPC. 



The simulations were done on matrix sizes 60 x 60 (DMPC- 
d54:DSPC 100:0, 70:30, 60:40, 50:50, 40:60, 30:70), 100 x 
100 (90:10, 80:20, 20:80, 10:90) and 350 x 350 (DSPC) ©. 
The system was equilibrated for 5000 or around the heat ca- 
pacity maxima for 10000 MC cycles. Simulations were con- 
ducted over 50000 or 100000 MC cycles, respectively. 



Results 

In this paper we discuss melting processes in binary lipid 
mixtures. Thereby, we focus on the melting of the two sin- 
gle components independently. Experimental FTIR and DSC 
measurements were accompanied by numerical simulations 
and an analysis of the corresponding phase diagram. 



MC simulations were used to evaluate a simple model based 
on the assumption of two lipid chain states, the considera- 
tion of nearest neighbor interactions only and the application 
of a hexagonal lattice. This model has already been shown 
to describe melting processes in DMPC DSPC mixtures suc- 
cessfully (I54ll57ll59h . The evaluation needs the determination 
of ten parameters. These are twice enthalpy and entropy dif- 
ferences between the two disordered and ordered chain states 
and six unlike nearest neighbor interaction parameters. The 
enthalpy and entropy changes can be deduced easily from 
heat capacity measurements of the single one component sys- 
tems. Two interaction parameters describe the cooperativ- 
ity, meaning the width of these transitions (I65[ |66J). The 



four parameters remaining need to be deduced from compar- 
isons with Cp-profiles obtained for different lipid composi- 
tions. Previously, we have determined the parameters for a 
DMPCDSPC system d59l) . In this study, however, we ex- 
changed the DMPC lipids with deuterated DMPC (DMPC- 
d54) which gives heat capacity profiles which are broadened 
and shifted by — 3.9°C to lower temperatures in compari- 
son to the nondeuterated DMPC (data not shown and see 
SB). This means that we had to determine the enthalpy 
and entropy changes and the cooperativity parameter of pure 
DMPC-d54. In principle using this model one needs to reestab- 
lish all parameters for any binary lipid system. In this work, 
however, we assumed that the unlike nearest neighbor inter- 
action parameters do not change by the replacing of DMPC 
with its deuterated analog and we adopted them from previ- 
ous simulations of the nondeuterated DMPCDSPC system 
rt59l) . To validate this, DSC measurements on an equimolar 
DMPCDSPC mixture were performed and they were com- 
pared to results simulated. We found that the values calcu- 
lated describe the melting profile measured well (see fig. [TJ. 
Therefore, we were convinced that it was feasible to further 
use the model with the interaction parameters given in tableQ] 

The further scope of this paper was to compare the temper- 
ature dependence of the simulated disordered chain ratios of 
DMPC-d54 and DSPC lipids with the experimental analysis 
of the melting processes using FTIR spectroscopy. This in- 
directly allowed us to compare the experimental results ob- 
tained with DSC and FTIR. 

In FTIR spectroscopy the absorption of infrared light due 
to vibrational degrees of freedom in the sample is studied. 
In lipid suspensions different vibrational modes of the lipids 
contribute to the absorption spectrum. Prominent examples 
used in our study are CH 2 or CD 2 (in the case of the deuter- 
ated DMPC lipids) stretch vibrations of the fatty acid lipid 
chains. The shape of the absorption band is temperature de- 
pendent as can be seen in the left panel of fig. [2] Example 
spectra of the symmetric and asymmetric CH 2 vibrational 
spectra are given for different temperatures of a mixture of 
DMPC-d54:DSPC 40:60. Deuterium atoms are heavier than 
hydrogen atoms. This leads to a shift of the CD 2 stretch vi- 
brations to lower wavenumbers and allows to distinguish the 
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two different lipid species in a binary lipid mixture. 



Comparison of FTIR Results and Monte Carlo Simula- 
tions 

In an investigation of structural changes in cytochrome c 
Heimburg and Marsh d67l) focused on temperature depen- 
dent changes of difference spectra. This was motivated by 
the idea that the corresponding amide I band is a convolution 
of a variety of different bands. The temperature dependence 
of each of the single bands differs and a complex tempera- 
ture dependent behavior is present. Therefore, we decided not 
only to investigate the temperature dependence of the absorp- 
tion maxima as done in the literature previously (28, 68, 69), 
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Figure 1 : Enthalpy and entropy changes and cooperativ- 
ity parameters of the single DMPC-d54 and DSPC lipid 
membranes were determined from heat capacity profiles. 
The remaining four unlike species nearest neighbor in- 
teraction parameters ofDMPC-d54:DSPC mixtures were 
assumed to equal the ones of a DMPCDSPC system 
$5$) . Simulating an equimolar mixture ofDMPC-d54 and 
DSPC gives a melting profile (markers) which describes 
well the measured heat capacity profile ( solid curve). Be- 
low three representative snapshots from simulations of 
the same mixture with a matrix size of 80 x 80 lipid chains 
are shown. The red and green colors represent ordered or 
disordered lipid chains, respectively. The lighter colors 
correspond to the DMPC-d54 lipid chains and the darker 
ones to the DSPC lipid chains. The snapshots were taken 
at three different temperatures. At a temperature of 15 
C the membrane is in the solid ordered phase while at a 
temperature of 55° C it is in the liquid disordered phase. 
In both cases lipids do not mix ideally, but lateral het- 
ereogeneities are present. Two macroscopic phases are 
present at a temperature o/33 C. 



but also to use a Difference Spectra method too. Further we 
studied the peak area of these difference spectra as a function 
of temperature. 

In the left panel of fig.[2]example spectra of the CH-2 stretch 
vibrations are given as a function of temperature of a DMPC- 
d54:DSPC 40:60 mixture. Normalizing the area of the spec- 
trum to 1 x cmT 1 and defining a reference spectrum which 
lies at a temperature below the melting transition regime, dif- 
ference spectra were calculated. Representative examples are 
displayed in the right panel of fig. |2] To follow the melt- 
ing transition the change in difference spectra peak area of 
the CH2 or the CD2 symmetric, asymmetric or both vibra- 
tional modes can be monitored. Our experiments showed that 
there is no difference between these three approaches (data 
not shown). 

The temperature dependent evolution of the difference spec- 
tra peak areas monitored are given in fig. [3] as solid curves 
for four different DMPC-d54:DSPC mixtures: (A) 70:30, 
(B) 50:50, (C) 40:60 and (D) 30:70. The panels on the left 
side depict the results related to the deuterated DMPC-d54 
lipids and the right ones give the data belonging to the DSPC 
lipids. The data obtained is also compared to simulation re- 
sults of the ratio of disordered lipid chains for each of the sin- 
gle species. In fig.[3]these are given as open circles (DMPC- 
d54) or squares (DSPC). In all instances we found an in- 
crease of the area with higher temperatures. In all cases a 
good agreement between the simulated and measured values 
is present. The deviations of the data simulated from the mea- 
sured ones are smaller in the cases of the DMPC-d54 lipids. 
This might be due to the fact that the determination of the 
model parameters needed is imperfect. More extended and 
time-consuming simulations might improve the data, too. 

As already mentioned above, previously melting processes 
of lipid membranes investigated by FTIR were analyzed mon- 
itoring the temperature dependence of the position of the ab- 
sorption maximum of either the symmetric or asymmetric 
stretch vibrations of CH2 and CD2. We analyzed the tem- 
perature dependence of the absorption maxima for several 
DMPC-d54:DSPC mixtures as earlier done by Leidy et al. 
J28I) . too. The results obtained are given by the dashed curves 
in fig. El 



We confirmed what has been reported by Leidy et al. (1281) 
previously. There are instances when the positions of the 
maxima of the symmetric stretch vibrations of the DSPC lipids 
do not increase with higher temperatures, but that for certain 
lipid mixture fractions they decrease, approach a minimum 
and then they start to become larger again. The minimum 
wavenumber is in all cases lower than the wavenumbers de- 
termined at lower temperatures. This behavior is in disagree- 
ment with the evolution of the difference spectra peak area 
and the MC simulation results. 

It is well known that macroscopic, but also local fluctuations 
of various membrane properties are enhanced in the vicinity 
of lipid membrane melting transitions (I56ll57l) . The determi- 
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Figure 2: 77;e /i? absorption spectrum contains different absorption maxima originating from varying vibrations. 
In the left panel the asymmetric and the symmetric CH2 stretch vibrations as measured in a sample of DMPC- 
d54:DSPC 40:60 at different temperatures are displayed. The absorption maxima of these are temperature de- 
pendent and lie around 2852 cmT 1 (symmetric) and at 2925 cmT x (asymmetric). The asymmetric and symmetric 
CD2 stretch vibrations are shifted to lower wavenumbers due to the heavier deuterium atoms (data not shown). 
In previous studies from other labs the melting processes were investigated following the temperature dependence 
of the absorption maxima. In this study, however, the melting process was also studied by defining a reference 
spectra and subtracting this one from all other spectra where in each cases the area was normalized to 1 x cm -1 
(Difference Spectra method). Example results of the difference spectra of the left panel are displayed in the right 
panel. As a measure of the progress of the melting process the area of one of the peaks was taken. 



nation of the derivatives of the data of the curves from fig. [3] 
as done by Leidy and collaborators) 28) provides a measure 
of the strength of fluctuations with respect to the single lipid 
components. In fig.|4]these derivatives are shown for the data 
simulated. The higher the value of the derivative the stronger 
the fluctuations. DMPC-d54 melts at lower temperatures than 
DSPC. This is why the corresponding derivatives in the left 
panels of fig.|4]have maxima at lower temperatures. This tem- 
perature, however, depends on the molar fraction of DSPC 
lipids and is higher with an increased DSPC lipid fraction. 
This is similar for the DSPC lipids only with the difference 
that the presence of DMPC-d54 lowers the temperatures at 
which fluctuations of the DSPC lipid chains are the strongest. 
In this context it should be noted that the component which is 
in the minority always displays two maxima. One of the max- 
ima is a local maximum at which fluctuations are enhanced, 
but they are lower than at the global maximum. The local 
maximum is found at a temperature at which the other com- 
ponent displays its strongest fluctuations. At an equimolar 
ratio of DMPC-d54:DSPC this is true for both components. 



The melting transition of the single lipid components happens 
over a broad temperature regime. 

The Phase Diagram 

Phase diagrams are generally constructed to easily access 
the phase separation behavior of complex lipid mixtures 
In our case a phase diagram depending on lipid composi- 
tion and temperature was constructed. A phenomenological 
method to construct these is in aligning tangents on the lower 
and upper temperature limits of heat capacity profiles. This 
results in the construction of a solidus and a liquidus line as 
indicated by open squares and the grey curves in fig. [5] If the 
lipid composition is known, the physical state of the mem- 
brane can be determined. The data obtained from the lower 
temperature limit defines the solidus line, whereas the one 
from the upper temperature limit gives the liquidus line. Be- 
low the solidus line the lipid membrane is in the solid ordered 
phase, while above the liquidus line it is in the liquid disor- 
dered phase. Therewith, the liquidus and the solidus lines 
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Figure 3: The FTIR results are compared with the Monte Carlo simulations. The data on four different DMPC- 
d54:DSPC ratios 70:30 (A), 50:50 (B), 40:60 ( C) and 30:70 (D) are displayed. In each panel the disordered chain 
ratios obtained from the simulations (open circles (DMPC-d54) and open squares (DSPC)) are compared with 
the changes of the peak area (solid curves) and the evolution of the position of the absorption maximum (dashed 
curves). In the latter case the antisymmetric CD2 (DMPC-d54) and the symmetric CH2 stretch vibrations were 
followed. The panels on the left side depict the melting process of the deuterated DMPC-d54 lipids and the right 
ones the DSPC lipid melting. The Monte Carlo simulations allow to combine both the FTIR and the calorimetric 
results. They provide a mean to relate the macroscopic with the microscopic behavior of the binary lipid membrane 
melting. 



define the phase coexistence regime too. The three snapshots 
in fig. Q] provide a visual impression of this. At temperatures 
of 15° C and 55° C the lipid membrane is in either the solid 
ordered or liquid disordered phase, respectively. It should, 
however, be mentioned that the mixing of the two lipids is 
in neither of the two cases ideal. The single lipid species 
cluster. At a temperature of 33° C the two phases coexist 
and a macroscopic phase separation is present. In an earlier 
publication we could show that discussing domain formation 
processes in lipid membranes one needs to distinguish micro- 
scopic and macroscopic phase separation and local fluctua- 
tions in chain state enter the physical picture. These fluctu- 
ations are enhanced at domain or phase boundaries. Small 
domains are also subject to strong fluctuations (|57I) . Analyz- 
ing phase diagrams one should note that these physical effects 
are in general not considered in its interpretation. Further, it 
has been pointed out that the construction of a phase diagram 
is strictly speaking only valid if the transition is of first-order 
J70I) . The transition of binary lipid mixtures of DMPC and 



DSPC has been claimed to be of second order and therefore 
the interpretation of the phase diagram has to be taken with 
care d55l) . 

With the determination of the derivatives of the curves in 
fig. [3] we were able to deduce the temperatures at which the 
fluctuations of the single components were enhanced. The 
lower temperature melting component, in this case DMPC- 
d54, has its maximum at lower temperatures than the one 
which melts at higher temperatures (DSPC). These tempera- 
tures, however, are higher or respectively lower than the melt- 
ing temperatures of the pure single components. The values 
are given in fig. [5] Open circles refer to the values originating 
from the simulations while the crosses come from the experi- 
ments. The temperature points determined using the two dif- 
ferent methods to analyze the FTIR absorption spectra are the 
same (data not shown). The experimental and the simulated 
data agree well. 

The simulations further allowed us to compare the tempera- 
tures at which the heat capacity curves are in a maximum and 
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Figure 4: The derivatives of the curves of the fraction of disordered DMPC-d54 and DSPC lipid chains as shown 
in fig.\3\are displayed in this figure. In detail these are the mixtures 70:30 (A), 50:50 (B), 40:60 ( C) and 30:70 (D). 
The derivatives belonging to DMPC-d54 are indicated by open circles, while the ones of DSPC by open squares. 
A maximum in the derivative means that at this temperature fluctuations are enhanced. 



the temperatures at which the derivatives of the disordered 
chain ratios display their maxima. We found that these tem- 
peratures equal each other {data not shown). In some cases 
only one heat capacity maximum is resolvable. The analysis 
of the single components, however, allows to distinguish the 
maxima of the melting of single components in all instances. 
The maxima of the component which is in a majority then 
agrees with the maximum of the heat capacity curve. 

In total, we showed a combined FTIR and MC simulation 
study. FTIR absorption spectra were analyzed with two meth- 
ods. One followed the temperature dependence of the absorp- 
tion maxima of the symmetric or asymmetric CH2 or CD2 
stretch vibrations. The other studied temperature dependent 
changes of the peak area of the difference spectra. Monitor- 
ing the positions of band maxima leads to incorrect results, 
but the simulated data describe well the curve shapes obtained 
with the Difference Spectra method. Our simulations describe 
well the evolution of the area changes. Melting of single 
components cannot be investigated in DSC. DSC alone only 
shows the bulk behavior. In FTIR spectroscopy it is possible 
to separate the lipid species individually. Therefore, using 
MC simulations we were able to link the two experimental 



approaches. 

Discussion 

In this paper we presented a combined numerical and experi- 
mental study exploring melting processes in binary lipid mix- 
tures. In detail, we focused on the melting of the single lipid 
components DMPC-d54 and DSPC independently. Monte 
Carlo simulations linked FTIR and DSC measurements to 
each other and therewith MC simulations proved another time 
to be a useful and powerful tool in exploring melting transi- 
tions in lipid membranes. 

FTIR, DSC and MC Simulations 

FTIR allowed us to measure the absorption spectra of lipid 
suspensions consisting of MLVs with varying DMPC-d54 and 
DSPC ratios in dependence of temperature. In particular we 
probed the CH2 and CD2 vibrational modes of the lipid fatty 
acid chains. The deuterium is heavier than the hydrogen atom 
which is why the symmetric and asymmetric bands of the 
CH 2 and CD 2 stretch vibrations display absorptions at dif- 
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ferent wavenumbers. This made it possible to distinguish the 
DMPC-d54 and DSPC lipids. 

We applied two different methods to analyze the absorption 
spectra measured. One was based on monitoring the temper- 
ature dependence of absorption maxima corresponding to the 
CH2 and CD2 stretch vibrations, respectively. The other one 
was built on following changes in the peak area of difference 
spectra (Difference Spectra method). We compared the cor- 
responding experimental results to the temperature dependent 
ratio of disordered chains of the single lipid species using MC 
simulations of a simple model (see fig. 01. 

These numerical simulations were based on a Doniach model 
of the lipid chain melting process assuming one ordered and 
one disordered lipid chain state. Nearest neighbor interac- 
tions only were considered and the lipid chains were arranged 
on a triangular lattice. Even though it is a minimalistic model 
which neglects that lipid chains may reflect different degrees 
of disorder and the non-consideration of a loss of lattice order 
during the melting process the calculated fraction of disor- 
dered lipid chains of the single lipid species describes the ex- 
perimental data obtained with the Difference Spectra method 




40 60 80 

fraction DSPC (%) 



100 



Figure 5: Phase diagrams allow to obtain an under- 
standing of the phase behavior of lipid mixtures in depen- 
dence of lipid ratios, temperature, pH or other outer pa- 
rameters. In this study binary lipid mixtures were used to 
study melting transitions in dependence of temperature. 
Therefore, a phase diagram with the fraction of DSPC 
lipids as the abscissa and temperature as the ordinate 
was constructed. The open squares and the grey lines 
represent the lower and upper temperature limits of the 
melting process as determined from heat capacity profiles 
simulated using the tangent method. Further the peak po- 
sitions of the derivatives of the disordered chain fractions 
simulated (open circles) and of the experimentally deter- 
mined area changes (crosses) were analyzed. The sim- 
ulated peaks determined from the derivatives agree with 
the maxima of the heat capacity profiles ( data not shown ) 



well. 

Previously this model has been shown to be suited to de- 
scribe heat capacity profiles as measured with DSC (154 157 , 
591) . Ten unique parameters are needed which were obtained 
from experimental heat capacity profiles. Six of them fol- 
low directly from the heat capacity curves of the one compo- 
nent systems. These are twice entropy and enthalpy changes 
and the corresponding cooperativity parameters. The missing 
four parameters have to be found by fitting the heat capacity 
profiles of the binary mixtures. They must be and they are 
the same for all of the possible lipid ratios. Due to the errors 
in the c p -profile measurements the parameter determination 
is subject to deviations from the ideal set of parameters. In 
our case we have taken the remaining four parameters from 
the DMPCDSPC system IH). This proved to be valid, but 
small deviations due to this assumption cannot be excluded, 
too. 

Heat capacity profiles contain macroscopic information on 
the melting behavior such as changes in enthalpy and entropy. 
However, DSC does not allow to determine the microscopic 
behavior of the melting process. Additionally, it does not al- 
low to investigate the melting of the single components. It 
is only suited to measure the bulk melting behavior. FTIR, 
however, provides a mean of studying the melting of single 
lipid species and observes the microscopic details, but lacks 
the possibility to probe macroscopic properties. The model 
underlying the MC simulations describes heat capacity pro- 
files and predicts the findings from FTIR well. Therefore, the 
MC simulations provide a mean to link the information of the 
two experimental techniques. 



The Melting Process 

The melting process of single components is adequately de- 
scribed using the numerical simulations. From the experi- 
mental and the simulated results and their corresponding deriva- 
tives it becomes clear that the lower temperature melting com- 
ponent (DMPC-d54) melts also in the binary mixture at a 
lower temperature than do the DSPC lipids (fig. [4j. Still, 
both lipid species melt over the whole temperature regime 
and not only at distinct temperatures at which fluctuations as- 
signed to the single lipid components are enhanced. In gen- 
eral, the component which is in the minority shows two max- 
ima. At one of these temperatures the strength of the fluctua- 
tions of the component which is in the majority is the highest. 
This effect underlines the cooperativity of the melting pro- 
cess. The melting of one lipid increases the probability of the 
melting of another lipid. 

The temperatures at which the fluctuations of the two com- 
ponents is strongest as obtained from the simulations and the 
experiments, agree well (fig. |5). Further, the numerical 
simulations allowed to calculate the heat capacity values and 
the disordered chain fractions simultaneously. Therefore, we 
were able to compare the temperatures at which the heat ca- 
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pacity is in its maximum with the temperatures determined 
from the derivatives of the disordered chain fractions. These 
temperatures agree with each other. In fig.[T]the heat capacity 
profile of the equimolar mixture displays two maxima. The 
maximum at lower temperature corresponds to the enhance- 
ment in fluctuations of the DMPC-d54 lipids and the one at 
higher temperatures to the one of the DSPC lipids. However, 
there are instances at which the heat capacity profile displays 
only one maximum even though one can determine two tem- 
peratures using the evolution of the single lipid disordered 
chain ratios. This is because DSC describes the bulk behav- 
ior and single events might not be detectable since they show 
only a small contribution to the overall melting. This needs 
attention in discussing melting processes in biological mem- 
branes. In the case of biological membranes it might even be 
difficult to distinguish the melting processes from the baseline 
of the calorimetric scan. This means that melting is present, 
but it cannot be detected at all using DSC. 

FTIR: The Two Approaches 

In this paper we used two different techniques to analyze 
the FTIR data. These were following the temperature depen- 
dence of the absorption maxima of either the symmetric CH2 
or the asymmetric CD2 stretch vibrations and the area change 
of one of the peaks of a difference spectra. 

Using the first method we found , in accordance to (28), that 



domains. Even if a membrane should be completely in the 
liquid disordered phase, lipids might not mix ideally and a 
heterogeneity exists as seen in fig.[T]and as already discussed 
elsewhere (1571) . This statement is also true for membranes in 
the liquid ordered phase (74J). The role of lipid membrane 
hetero gene ities has been discussed exhaustingly in the liter- 



the absorption maximum of the CH2 stretch vibrations of the 
DSPC lipids developed towards smaller wavenumbers when 
the DMPC-d54 started to display stronger fluctuations (see 
the right panels of fig- [3j. It has been interpreted that in these 
cases the DSPC lipids start to order. Concluding from the 
data presented in this work and by Leidy et al. (1281) this or- 
dering should be stronger than in the presence of a complete 
solid ordered membrane. The wavenumber at which the ab- 
sorption is at a maximum drops to wavenumbers even smaller 
than at low temperatures. There is, however, no further sup- 
porting evidence that lipids might be able to order stronger 
than at low temperatures. This behavior is in disaccord of 
what we have found with the second method and the Monte 
Carlo simulations. Therefore, it seems reasonable that the 
Difference Spectra method yields results which describe the 
melting process correctly and monitoring band shifts only can 
lead to incorrect results. 



Biological Relevance 

The model of a biological membrane by Singer and Nicol- 
son considers the biological membrane to be homogeneous 
(fnl) . In the recent years this view has changed. A hetero- 
geneous membrane model has overtaken this older picture 
(172l l73l) . Different kinds of heterogeneities can be present 
in biological membranes, such as protein aggregates or lipid 



ature Q |M [H, S3 EH Q SslS). Thereby, 



it has not only been pointed out that the coexistence of solid 
ordered and liquid disordered phases might be important (H3, 
13 41 , 43 , 44j), but also the enrichment of single lipi d species 
in certain regions of the lipid membrane d42l I45LI461 \4m . 

Melting of single components cannot be detected from the 
bulk melting in the calorimeter. However, most of the pub- 
lications cited here deduce phase separation of natural mem- 
branes from calorimetric results. Therefore, details of these 
processes might not be revealed. It is also possible that the 
melting process in a biological membrane displays a broad 
profile and the heat changes detected by a calorimeter cannot 
be distinguished from the background signal. 

In this context not only the various lipid membrane phases 
have to be named, but also the possibility of enrichment of 
certain lipid species in the certain lipid domains. Therefore, 
however, a detailed understanding of lipid melting transitions 
and lateral structuring need to be obtained. The existence of 
heterogeneities depends on the thermodynamical behavior of 
the different lipid species. As seen in this study they melt over 
a broad temperature regime with particular temperatures at 
which they display strong fluctuations in more complex lipid 
mixtures. These temperatures depend on the lipid species, 
but also on its interactions with other kinds of lipid present. 
Melting processes might not be finished even though a heat 
capacity profile might not indicate any such events. The sen- 
sitivity might be not sufficient. FTIR spectroscopy, however, 
provides means to distinguish the melting of single compo- 
nents. 

A possible reason for the big variety in lipid species might be 
to ensure a needed lipid membrane heterogeneity and there- 
with to play a role membrane function. This then provides an 
explanation why a dependence of lipid synthesize in biologi- 
cal membranes on outer physical parameters such as pressure 
or temperature is required. 
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